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ABSTRACT 

The Planck design and scanning strategy provide many levels of redundancy that can be exploited to provide tests of internal consistency. One of 
the most important is the comparison of the 70 GHz (amplifier) and 100 GHz (bolometer) channels. Based on different instrument technologies, 
with feeds located differently in the focal plane, analysed independently by different teams using different software, and near the minimum of 
diffuse foreground emission, these channels are in effect two different experiments. The 143 GHz channel has the lowest noise level on Planck, and 
is near the minimum of unresolved foreground emission. In this paper, we analyse the level of consistency achieved in the 2013 Planck data. We 
concentrate on comparisons between the 70, 100, and 143 GHz channel maps and power spectra, particularly over the angular scales of the first 
and second acoustic peaks, on maps masked for diffuse Galactic emission and for strong unresolved sources. Difference maps covering angular 
scales from 8° to 15' are consistent with noise, and show no evidence of cosmic microwave background structure. Including small but important 
corrections for unresolved-source residuals, we demonstrate agreement (measured by deviation of the ratio from unity) between 70 and 100 GHz 
power spectra averaged over 70 < ^ < 390 at the 0.8 % level, and agreement between 143 and 100 GHz power spectra of 0.4 % over the same I 
range. These values are within and consistent with the overall uncertainties in calibration given in the Planck 2013 results. We also present results 
based on the 2013 likelihood analysis showing consistency at the 0.35 % between the 100,143, and 217 GHz power spectra. We analyse calibration 
procedures and beams to determine what fraction of these differences can be accounted for by known approximations or systematic errors that 
could be controlled even better in the future, reducing uncertainties still further. Several possible small improvements are described. Subsequent 
analysis of the beams quantifies the importance of asymmetry in the near sidelobes, which was not fully accounted for initially, affecting the 
70/100 ratio. Correcting for this, the 70, 100, and 143 GHz power spectra agree to 0.4 % over the first two acoustic peaks. The likelihood analysis 
that produced the 2013 cosmological parameters incorporated uncertainties larger than this. We show explicitly that correction of the missing near 
sidelobe power in the HRI channels would result in shifts in the posterior distributions of parameters of less than 0.3cr except for As, the amplitude 
of the primordial curvature perturbations at 0.05 Mpc"E which changes by about Icr. We extend these comparisons to include the sky maps from 
the complete nine-year mission of the Wilkinson Microwave Anisotropy Probe (WMAP), and find a roughly 2 % difference between the Planck and 
WMAP power spectra in the region of the first acoustic peak. 

Key words. Cosmology: observations — Cosmic background radiation — Instrumentation: detectors 


1. Introduction 


* Corresponding author: C. R. Lawrence Charles. lawrence@jpl. 
nasa.gov 


This paper, one of a set associated with the 2013 release of data 
from the Planck^ mission (Planck Collaboration I 2014), de- 

^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
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scribes aspects of the internal consistency of the Planck data in 
the 2013 release not addressed in the other papers. The Planck 
design and scanning strategy provide many levels of redundancy, 
which can be exploited to provide tests of consistency (Planck 
Collaboration I 2014), most of which are carried out routinely 
in the Planck data processing pipelines (Planck Collaboration II 
2014; Planck Collaboration VI 2014; Planck Collaboration XV 
2014; and Planck Collaboration XVI2014). One of the most im¬ 
portant consistency tests for Planck is the comparison of the 
LFI and HFI channels, and indeed this was a key feature of 
its original experimental concept. Based on different instrument 
technologies, with feeds located differently in the focal plane, 
and analysed independently by different teams, these two instru¬ 
ments provide a powerful mutual assessment and test of system¬ 
atic errors. This paper focuses on comparison of the LFI and HFI 
channels closest in frequency to each other and to the diffuse 
foreground minimum, namely the 70 GHz (LFI) and 100 GHz 
(HFI) channels^, together with the 143 GHz HFI channel, which 
has the greatest sensitivity to the cosmic microwave background 
(CMB) of all the Planck channels. 

Quantitative comparisons involving different frequencies 
must take into account the effects of frequency-dependent fore¬ 
grounds, both diffuse and unresolved. Planck processing for the 
2013 results proceeds along two main lines, depending on the 
scientific purpose. For non-Gaussianity and higher-order statis¬ 
tics, and for the ^ < 50 likelihood, diffuse foregrounds are sep¬ 
arated at map level (Planck Collaboration XII 2014). Only the 
strongest unresolved sources, however, can be identified and 
masked from the maps, and the effects of residual unresolved 
foregrounds must be dealt with statistically. They therefore re¬ 
quire corrections later in processing (e.g., Planck Collaboration 
XXIII2014, Planck Collaboration XXIV 2014). For power spec¬ 
tra, the ^ > 50 likelihood, and parameters, both diffuse and un¬ 
resolved source residuals are handled in the power spectra with 
a combination of masking and fitting of a parametric foreground 
model (Planck Collaboration XVI 2014). 

In this paper, we compare Planck channels for consistency 
in two different ways. First, in Sects. 2 and 3, we compare fre¬ 
quency maps from the Planck 2013 data release - available from 
the Planck Legacy Archive (PLA)^ and referred to hereafter as 
PLA maps - and power spectra calculated from them, looking 
first at the effects of noise and foregrounds (both diffuse and un¬ 
resolved), and then, in Sect. 4, at calibration and beam effects. 
This comparison based on publicly released maps ties effects in 
the data directly to characteristics of the instruments and their 
determination. This examination has provided important insights 
into our calibration and beam determination procedures even 
since the 2013 results were first released publicly in March 2013, 
confirming the validity of the 2013 cosmological results, resolv¬ 
ing some issues that had been contributing to the uncertainties, 
and suggesting future improvements that will reduce uncertain¬ 
ties further. 

Second, in Sect. 5, we compare power spectra again, this 
time from the “detector set” data at 100, 143, and 217 GHz (see 
Table 1 in Planck Collaboration XV 2014) used in the likelihood 


entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci¬ 
entific consortium led and funded by Denmark. 

^ The frequency at which extragalactic foregrounds are at a minimum 
level depends on angular scale, shifting from around 65 GHz at low I to 
143 GHz at ^ ^ 200. 

^ http://archives.esac.esa.int/pla2 


analysis described in Planck Collaboration XV (2014), but ex¬ 
tending that analysis to include 70 GHz as well. These detector- 
set/likelihood comparisons give a measure of the agreement be¬ 
tween frequencies in the data used to generate the Planck 2013 
cosmological parameter results in Planck Collaboration XVI 
(2014). Taking into account differences in the data and process¬ 
ing, the same level of consistency is seen as in the comparison 
based on PLA frequency maps in Sect. 3. We then show that the 
small changes in beam window functions discussed in Sect. 4 
have no significant effect on the 2013 parameter results other 
than the overall amplitude of the primordial curvature perturba¬ 
tions at 0.05 Mpc“\ As. 

After having established consistency within the Planck data, 
specifically agreement between 70, 100, and 143 GHz over the 
first acoustic peak to better than 0.5 % in the power spectrum, 
in Sect. 6 we compare Planck with WMAP, specifically the 
WMAP9 release^. The absolute calibration of the Planck 2013 
results is based on the “solar dipole” (i.e., the motion of the 
Solar System barycentre with respect to the CMB) determined 
by WMAPl (Hinshaw et al. 2009), whose uncertainty leads to a 
calibration error of 0.25 % (Planck Collaboration V 2014). For 
the Planck channels considered in this paper, the overall cali¬ 
bration uncertainty is 0.6% in the 70 GHz maps and 0.5 % in 
the 100 and 143 GHz maps (1.2% and 1.0%, respectively, in 
the power spectra; Planck Collaboration I 2014, Table 6). When 
comparing Planck and WMAP calibrated maps, however, one 
should remove from these uncertainties in the Planck maps the 
0.25 % contribution from the WMAP dipole, since it was the ref¬ 
erence calibrator for both LFI and HFI. In the planned 2014 re¬ 
lease, the Planck absolute calibration will be based on the “or¬ 
bital dipole” (i.e., the modulation of the solar dipole due to the 
Earth’s orbital motion around the Sun), bypassing uncertainties 
in the solar dipole. 

Throughout this paper we refer to frequency bands by their 
nominal designations of 30, 44, 70, 100, 143, 217, 353, 545, and 
857 GHz for Planck and 23, 33, 41,61, and 94 GHz for WMAP\ 
however, we take bandpasses into account in all calculations. 
The actual weighted central frequencies determined by convo¬ 
lution of the bandpass response with a CMB spectrum are 28.4, 
44.1, 70.4, 100.0, 143.0, 217.0, 353.0, 545.0, and 857.0 GHz for 
Planck, and 22.8, 33.2, 41.0, 61.4, and 94.0 GHz for WMAP. 
These correspond to the effective frequencies for CMB emis¬ 
sion. For emission with different spectra, the effective frequency 
is slightly shifted. 

The maps discussed in this paper are structured accord¬ 
ing to the HEALPix^ scheme (Gorski et al. 2005) displayed in 
Mollweide projections in Galactic coordinates. 


2, Comparison of frequency maps 

The Planck 2013 data release includes maps based on 
15.5 months of data, as well as maps of subsets of the data that 
enable tests of data quality and systematic errors. Examples in¬ 
clude (see Planck Collaboration I 2014 for complete descrip¬ 
tions) single survey maps and half-ring difference maps, made 
by splitting the data from each pointing period of the satellite 
into halves, making separate sky maps from the two halves, and 
taking the difference of the two maps. Half-ring maps are partic¬ 
ularly useful in characterizing the noise, and also enable signal 


Available from the LAMBDA site: http: //lambda. gsfc. nasa. 

gov 

^ See http://healpix.jpl■nasa.gov 
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Fig. 1. Sky masks used for spectral analysis of the Planck 
70, 100, and 143 GHz maps. The light blue, yellow, and red 
masks leave observable sky fractions /sky of 39.7 %, 59.6 %, 
and 69.4 %, respectively, and are named GAL040, GAL060, and 
GAL070 in the PLA. These masks are extended by exclusion of 
unresolved sources in the PCCS 70, 100, and 143 GHz source 
lists above the 5 cr flux density cuts. 


estimation based on cross-spectra, with signiflcant noise reduc¬ 
tion compared to auto-spectra. 

The 100 and 143 GHz maps are released at HEALPix reso¬ 
lution Aside = 2048, with Apix = 12 x 5 x 10^ pixels 

of approximately 1.'7. Although the LFI maps are generally re¬ 
leased at Aside = 1024, the 70 GHz maps are also released at 
A^side = 2048. All mapmaking steps except map binning at the 
given pixel resolution are the same for the two resolutions. In 
this paper we use the 70 GHz maps made at Aside = 2048 for 
comparison with the 100 and 143 GHz sky maps. 

2.1. Sky masks 

Comparison of maps at different frequencies over the full sky 
is quite revealing of foregrounds, as will be seen. For most 
purposes in this paper we need to mask regions of strong 
foreground emission. We do this using the publicly-released^ 
Galactic masks GAL040, GAL060, and GAL070, shown in 
Fig. 1. These leave unmasked /sky = 39.7 %, 59.6 %, and 69.4 % 
of the full sky. We mask unresolved (“point”) sources detected 
above 5cr in the 70, 100, and 143 GHz channels, as described 
in the Planck Catalogue of Compact Sources (PCCS; Planck 
Collaboration XXVIII 2014). The point source masks are cir¬ 
cular holes centred on detected sources with diameter 2.25 times 
the FWHM beamsize of the frequency channel in question. The 
masks are unapodized, as the efl'ect of apodization on large an¬ 
gular scales is primarily to improve the accuracy of covariance 
matrices. 


the few microkelvin level in the difference maps smoothed to 8° 
shown in Appendix A. In those cases, a small offset adjustment, 
typically no more than a few microkelvin, is made to keep the 
mean value very close to zero in patches of sky visually clear of 
foregrounds. 

2.3. Comparisons 

Figure 2 shows the monopole- and dipole-removed maps at 70, 
100, and 143 GHz, along with the corresponding half-ring dif¬ 
ference maps. Figure 3 shows the difference maps between these 
three frequencies. The strong frequency-dependence of fore¬ 
grounds is obvious. Equally obvious, and the essential point 
of the comparison, is the nearly complete nulling of the CMB 
anisotropies. This shows that these three channels on Planck are 
measuring the same CMB sky. 

For a quantitative comparison, we calculate root mean square 
(rms) values of unmasked regions of the frequency and dif¬ 
ference maps shown in Figs. 2 and 3, for the three masks 
shown in Fig. 1. To avoid spurious values caused on small 
scales by the differing angular resolution of the three frequen¬ 
cies, and on large scales by diffuse foregrounds, we first smooth 
the maps to a common resolution of 15'. We then smooth 
them further to 8° resolution, and subtract the 8° maps from 
the 15' maps. This leaves maps that can be directly compared 
for structure on angular scales from 8° to 15'. We calculate 
rms values for half-ring sum (“Freq.”) and half-ring difference 
(“Diff.”) maps at 70, 100, and 143 GHz, and for the frequency- 
difference maps 70 GHz-100 GHz, 70 GHz - 143 GHz, and 
100 GHz - 143 GHz. The rms values are given in Table 1. The 
maps are shown in Fig. A.l. Histograms of the frequency maps 
and difference maps are shown in Fig. 4. 

Except for obvious foreground structures and noise, the dif¬ 
ference maps lie close to zero, showing the excellent agreement 
between the three Planck frequencies for the CMB anisotropies. 

The map comparisons give a comprehensive view of consis¬ 
tency between 70, 100, and 143 GHz, but the two-dimensional 
nature of the comparisons makes it somewhat difficult to grasp 
the key similarities. To make this easier, we turn now to compar¬ 
isons at the power spectrum level. 


3. Comparison of power spectra from 2013 results 
frequency maps 

Power spectra of the unmasked regions of the maps are estimated 
as follows. 


2.2. Monopole/dipole removal 

The Planck data have an undetermined absolute zero level, and 
the Planck maps contain low-amplitude offsets generated in the 
process of mapmaking, as well as small residual dipoles that re¬ 
main after removal of the kinematic dipole anisotropy. We re¬ 
move the ^ = 0 and £ = \ modes from the maps using x^- 
minimization and the GAL040 mask, extended where applicable 
to a constant latitude of +45°. Diffuse Galactic emission at both 
low and high frequencies is still present even at high latitude, 
so this first step can leave residual offsets that become visible at 

^ Available from the Planck Legacy Archive: http://archives. 
esac.esa.int/pla2 


- Starting from half-ring maps (section 5.1 of Planck 
Collaboration I 2014), cross-spectra are computed on the 
masked, incomplete sky using the HEALPix routine anafast 
with Z’ = 0, 1 removal. These are so-called “pseudo-spectra.” 

- The MASTER spectral coupling kernel (Hivon et al. 2002), 
which describes spectral mode coupling on an incomplete 
sky, is calculated based on the mask used. The pseudo¬ 
spectra from the previous step are converted to 47r-equivalent 
amplitude using the inverse of the MASTER kernel. 

- Beam and pixel smoothing effects are removed from the 
spectra by dividing out the appropriate beam and pixel win¬ 
dow functions. Beam response functions in £ space are re¬ 
quired. We use the effective beam window functions derived 
using FEBeCoP (Mitra et al. 2011). 
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Fig. 2. Sky maps used in the analysis of Planck data consistency. Top row: 70 GHz. Middle row: 100 GHz. Bottom row: 143 GHz. 
Left column: signal maps. Right column: noise maps derived from half-ring differences. All maps are Aside = 2048. These are the 
publicly-released maps corrected for monopole and dipole terms as described in the text. The impression of overall colour differences 
between the maps is due to the interaction between noise, the colour scale, and display resolution. For example, the larger positive 
and negative swings between pixels in the 70 GHz noise map pick up darker reds and blues farther from zero. Smaller swings around 
zero in the 100 and 143 GHz noise maps result in pastel yellows and blues in adjacent pixels, which when displayed at less than 
full-pixel resolution give an overall impression of green, a colour not used in the colour bar. 


3.1. Spectral analysis of signals and noise 

Figure 5 shows the signal (half-ring map cross-spectra), and 
noise (half-ring difference map auto-spectra) of the 70, 100, 
and 143 GHz channels. As stated earlier, the 70-100 GHz chan¬ 
nel comparison quantifies the cross-instrument consistency of 
Planck. 

This description of the statistics of noise contributions to 
the empirical cross-spectra derived from the Planck sky maps 
sets up the analysis of inter-frequency consistency of Planck 
data. The pure instrumental noise contribution to the empiri¬ 
cal cross-spectra is very small over a large ^-range for the HFI 
channels, and at 70 GHz over the Arange of the first peak in 
the spectrum, where we now focus our analysis. Cosmic vari¬ 


ance is irrelevant for our discussion because we are assessing 
inter-frequency data consistency, and the instruments observe 
the same CMB anisotropy. Any possible departures from com¬ 
plete consistency of the measurements must be accounted for by 
frequency-dependent foreground emission, accurate accounting 
of systematic effects, or (at a very low level) residual noise. 

3.2. Spectral consistency 

Figure 6 shows spectra of the 70, 100, and 143 GHz maps for 
the three sky masks, differences with respect to the Planck 2013 
best-fit model, and ratios of different frequencies. In the ^-range 
of the first peak and below, the 143/100 ratio shows the ef- 
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Fig. 3. Difference maps. Top: 100 GHz minus 70 GHz. Middle: 143 GHz minus 100 GHz. Bottom: 143 GHz minus 70 GHz. All sky 
maps are smoothed to angular resolution FWHM = 15' by a filter that accounts for the difference between the effective beam 
response at each frequency and a Gaussian of FWHM 15'. These maps illustrate clearly the difference in the noise level of the 
individual maps, excellent overall nulling of the CMB anisotropy signal, and frequency-dependent foregrounds. The 100-70 
difference shows predominantly CO (7 = 0^1) emission (positive) and free-free emission (negative). The 143 - 100 difference 
shows dust emission (positive) and CO emission (negative). The 143 - 70 difference shows dust emission (positive) and free-free 
emission (negative). The darker stripe in the top and bottom maps is due to reduced integration time in the 70 GHz channel in the 
first days of observation (see Planck Collaboration II 2014, Sect. 9.5). 
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Table 1. Rms values of the unmasked regions of the frequency and difference maps shown in Fig. A.l, and for the /sky = 69.4 %, 
59.6%, and 39.7% masks shown in Fig. 1. Diagonal blocks give the rms values for half-ring sums (“Freq.”) and differences 
(“Diff.”) of single-frequency maps. Off-diagonal blocks give the same quantities for froquoncy-difference maps (70 GHz - 100 GHz, 
70 GHz - 143 GHz, and 100 GHz - 143 GHz. As described in the text, the maps are smoothed to a common resolution of 15', some¬ 
what lower than the resolution of the 70 GHz maps. In addition, structure on scales larger than 8 ° is determined and removed from 
all maps to avoid introducing biases from residual monopoles and dipoles, so that only structure from 15' to 8 ° in angular scale is 
included in these calculations. 


RMS [//K] 

70 GHz 100 GHz 143 GHz 


/sky Freq. Diff. Freq. Diff. Freq. Diff. 


39.7%. 90.44 28.62 29.01 28.93 29.00 28.69 

70 GHz 59.6%. 90.09 29.48 29.66 29.47 29.77 29.22 

69.4%. 90.12 29.46 29.79 29.36 30.03 29.12 

39.7%. 85.63 4.27 5.49 4.76 

100 GHz 59.6%. 85.05 4.38 6.09 4.83 

69.4%. 85.16 4.39 6.76 4.81 

39.7%. 85.70 2.11 

143 GHz 59.6%. 85.23 2.17 

69.4%. 85.45 2.18 


fects of residual diffuse foreground emission outside the masks. 
The largest mask reduces the detected amplitude, but does not 
remove it completely. The frequency dependence of the ratios 
conforms to what is well known, namely, that diffuse fore¬ 
ground emission is at a minimum between 70 and 100 GHz. The 
143/100 pair is more affected by diffuse foregrounds than the 
70/100 pair, as the dust emission gets brighter at 143 GHz. The 
effects of residual unresolved foregrounds in Fig. 6 are discussed 
in the next section. 

Near the first acoustic peak, measurements in the three 
Planck channels agree to better than one percent of the CMB 
signal, and to much better than their uncertainties, which are 
dominated by the effects of cosmic/sample variance (see Fig. 6 ). 
Inclusion of cosmic/sample variance is essential for making 
inferences about the underlying statistical processes of the 
Universe; however, since the receivers at all frequencies are 
observing a single realization of the CMB, cosmic variance is 
irrelevant in the comparison of the measurements themselves. 
Figure 7 is the same as the top two middle panels of Fig. 6 (i.e., 
over 60 % of the sky), but without inclusion of cosmic/sample 
variance in the uncertainties. As can be seen, cosmic/sample 
variance completely dominates the measurement uncertainties 
up to multipoles of 400, after which noise dominates. 

3.3. Residual unresolved sources 

Figure 6 shows that while diffuse foregrounds are significant for 
low multipoles, they are much less important on smaller angular 
scales. To see clearly the intrinsic consistency between frequen¬ 
cies, however, we must remove the effects of unresolved sources. 

Discrete extragalactic foregrounds comprise synchrotron ra¬ 
dio sources, Sunyaev-Zeldovich (SZ) emission in clusters, and 
dust emission in galaxies. These have complicated behaviour in 
/ and y. All have a Poisson part, but the SZ and cosmic infrared 
background sources also have a correlated part. These are the 
dominant foregrounds (for a 39.7 % Galactic mask) for € ^ 200. 
For frequencies in the range 70-143 GHz and multipoles in the 
range 50-200, they stay below 0.2 %. The minimum in unre¬ 
solved foregrounds remains at 143 GHz, with less than 2 % con¬ 
tamination up to ^ = 1000 . 


Discrete sources detected above 5 cr in the PCCS (Planck 
Collaboration XXVIII 2014) are individually masked, as de¬ 
scribed in Sect. 2. Corrections for residual unresolved radio 
sources are determined by fitting the differential Euclidean- 
normalized number counts jdS in Jy^'^ sr“^at each fre¬ 

quency with a double power law plus Euclidean term: 

^ 4 . 55/2 

5 ! dNIdS = [(5/5+^^(1 _g-s/5E)’ 

where Af is the amplitude at faint flux density levels, S 1 is the 
first faint flux density level, bn is the exponent of the first power 
law at faint flux densities, 5 2 is the second faint flux density 
level, bf 2 is the exponent of the second power law at faint flux 
densities, Ae is the amplitude of the Euclidean part, i.e., at large 
flux density, and ^e is the flux density level for the Euclidean 
part (^ 1 Jy). These are then integrated from a cutoff flux density 
corresponding to the 5 cr selection limit in the PCCS at 143 GHz, 
and the equivalent levels for a radio source with S oc y at 
100 and 70 GHz. Thermal SZ and GIB fluctuations are fitted 
as part of likelihood function determination described in Planck 
Collaboration XV (2014); the values found there are used here. 
Figure 8 shows the level of these corrections, while Fig. 9 shows 
the ratios of power spectra after the corrections are made. 

3.4. Assessment 

The 70/100 and 143/100 ratios in Fig. 9, for 59.6 % of the sky, 
averaged over the range 70 < ^ < 390 where the 70 GHz signal- 
to-noise ratio (S/N) is high, are 1.0080 and 1.0045, respectively. 
Over the range 70 < £ < 830, the ratios are 1.0094 and 1.0043, 
respectively. Table 2 collects these ratios and following ones for 
easy comparison. 

Section 7.4 of Planck Collaboration VI (2014) uses the 
SMICA code to intercalibrate on the common CMB anisotropies 
themselves, with results given in figure 35 of that paper. For 40 % 
of the sky, the 70/100 and 143/100 power ratios are 1.006 and 
1.002 over the range 50 < i < 300, and 1.0075 and 1.002 over 
the range 300 < i < 100. (These gain ratios from figure 35 
of Planck Collaboration VI 2014 must be squared for compari¬ 
son with the power ratios discussed in this section and given in 
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Fig. 4. Signal and noise for the frequency maps of Fig. 2 (left panel) and the difference maps of Fig. 3 (right panel), with the 59.6 % 
mask in all cases. The broader, signal+noise curves are nearly Gaussian due to the dominant CMB anisotropies. The 70 GHz curve 
is broader than the 100 and 143 GHz curves because of the higher noise level, but is still signal-dominated for \dTIT\ > 50pK. 
The narrower noise curves, derived from the half-ring difference maps, are not Gaussian because of the scanning-induced spatial 
dependence of pixel noise in Planck maps. The considerably higher noise level of the 70 GHz map is again apparent. The histograms 
of the difference maps show noise domination near the peak of each pair of curves (the signal-rnoise and noise curves overlap^ The 
pairs involving 70 GHz are wider and dominated by the 70 GHz noise, but the wings at low pixel counts show the signature of 
foregrounds that exceed the noise levels, primarily dust and CO emission in the negative wing, and free-free and synchrotron 
emission in the positive wing. In the low-noise 100 minus 143 GHz pairs, the signal, due mostly to dust emission in the negative 
wing and to free-free and CO residuals in the positive wing, stands out clearly from the noise. 


Table 2.) The SMICA equivalent power ratios are systematically 
about 0.2 % closer to unity than those calculated in this section; 
however, in broad terms the two methods give remarkably sim¬ 
ilar results. Moreover, the absolute gain calibration uncertain¬ 
ties given in Planck Collaboration V (2014, Table 8) and Planck 
Collaboration VIII (2014) are 0.62 % for 70 GHz and 0.54 % for 
100 GHz and 143 GHz. The agreement at the power spectrum 
level between 70, 100, and 143 GHz is quite reasonable in terms 
of these overall uncertainties. We will return to comparisons of 
spectra in Sects. 4 and 5. 

We are working continuously to refine our understanding of 
the instrument characteristics, implement more accurate calibra¬ 
tion procedures, and understand and control systematic effects 
better. All of these will lead to reduced errors and uncertainties 
in 2014. In the next section we describe an analysis of beams 
and calibration procedures that has already been beneficial. 

4. Beams, beam transfer functions, and calibration 

The residual differences that we see in Sect. 3 are small but not 
negligible. We now address the question of whether they may be 
due to beam or calibration errors. Detailed descriptions and anal¬ 
yses of the LFI and HFI beams and calibration are contained in 
Planck Collaboration IV (2014), Planck Collaboration V (2014), 
Planck Collaboration VII (2014), and Planck Collaboration VIII 
(2014). In this section we summarize our present understanding 


of calibration and beam effects for the two instruments, explain 
the reasons for the approximations that have been made in data 
processing, provide estimates for the impact of these approxima¬ 
tions on the resulting maps and power spectra, and outline plans 
for changes to be implemented in the 2014 data release. We will 
show that the small differences between LFI and HFI at interme¬ 
diate i seen in Fig. 9 are significantly reduced by improvements 
in our understanding of the near sidelobes in HFI, which affect 
the window functions in this i range. 


4.1. Beann definitions 

Calibration of the CMB channels (30 to 353 GHz) is based on the 
dipole anisotropies produced by the motion of the Sun relative to 
the CMB and of the modulation of this dipole by the motion of 
the spacecraft relative to the Sun (which we refer to as the solar 
and orbital dipoles, respectively). For the 2013 data release and 
all LFI and HFI frequency channels considered in this paper, the 
time-ordered data have been fit to the solar dipole as measured 
by WMAP (Hinshaw et al. 2009). The present analysis aims to 
show that the LFI-HFI differences at intermediate £ seen in Fig. 9 
are understood within the present uncertainties due to beams, 
calibration, and detector noise. Planck Collaboration IV (2014) 
and Planck Collaboration VII (2014) define three regions of the 


7 























Planck Collaboration: Planck 2013 results. XXXL 



Fig. 5. Planck 70, 100, and 143 GHz CMB anisotropy power spectra computed for the GAL060 mask. Mask- and beam-deconvolved 
cross-spectra of the half-ring maps show the signal; auto-correlation spectra of the half-ring difference maps show the noise. Points 
show single multipoles up to ^ = 1200 for 70 GHz and I = 1700 for 100 and 143 GHz. Heavy solid lines show = 20 boxcar 
averages. The S/N near the first peak {i = 220) is approximately 80, 1900, and 6000 for 70, 100, and 143 GHz, respectively. Noise 

power is calculated according to the large-^ approximation, i.e., as distribution with mean Q and rms Q [/sky(2^ + l)/2] ^ . 


Pairs of thin lines mark +3 cr bands of noise power around the noise spectra. We translate this statistical spread of noise power Qs 
into the signal spectra estimated via half-ring map cross-spectra. Under the simplifying assumption that each Q of the noise in 
the cross-spectrum at high-^ is distributed as a sum of {2i -r 1) products of independent Gaussian deviates, each with variance 
derived from the half-ring difference maps, the Gaussianized high-^ noise in the cross-spectra has zero mean and rms of 

1/2 


2 qoise 

2qoise 


[/sky (2^ +1)] • Pairs of thin lines mark +1 cr bands of noise around the boxcar-averaged cross-spectra. 


beam response (see Fig. 10, Fig. 1 of Planck Collaboration IV 
2014, and Fig. 5 of Tauber et al. 2010), as follows. 

The nominal beam or main beam is that portion used to create 
the beam window functions for the 2013 data release. The nom¬ 
inal beam carries most of the beam shape information and more 
than 99 % of the total solid angle, and therefore has most of the 
information needed for the 2013 cosmological analysis. The an¬ 
gle from the beam centre to the boundary of the nominal beam 
varies with frequency and instrument, and is 1?9, 1?3, and 0?9 
for 30, 44, and 70 GHz, respectively, and 0?5 for 100 GHz and 
above. 

The near sidelobes comprise any effective solid angle within 5° 
of the centre of the beam that is not included in the nominal 
beam. The response to the dipole from this region of the beam is 
very similar to that from the nominal beam, and unaccounted-for 
near sidelobe response leads to errors in the window function. 

The far sidelobes comprise the beam response more than 5° from 
the centre. Because of the geometry of the telescope and baffles, 
the bulk of this solid angle is at large angles from the line of 
sight, not far from the spin axis, and not in phase with the dipole 
seen in the nominal beam, and therefore has little effect on the 


dipole calibration. However, the secondary mirror spillover, con¬ 
taining typically 1/3 of the total power in the far sidelobes, is in 
phase with the dipole, and affects the calibration signal. The in¬ 
accuracy introduced by approximating the optical response with 
the nominal beam normalized to unity is corrected to first order 
by our use of a pencil (d-function) beam to estimate the calibra¬ 
tion. 

For reference, for the 2013 release we estimated a contribu¬ 
tion to the solid angle from near sidelobes of 0.08 %, 0.2 %, and 
0.2 %, and from far sidelobes of 0.62 %, 0.33 %, and 0.31 %, for 
70, 100, and 143 GHz, respectively, of which 0.12%, 0.075 %, 
and 0.055 % is from the secondary spillover referred to above. 
Recent analysis, detailed in Appendix C, has resulted in a new 
estimate for the near sidelobe contribution for 100 and 143 GHz 
of 0.30 + 0.2 % and 0.35 + 0.1 %, respectively^. The impact of 
this is described below. 


^ The solid angle statistical errors are 0.53 % and 0.14 % at 100 and 
143 GHz, respectively. 




Power ratio AC^ 1 )/27 t [/i-K^] C^^(^+1)/27 t [yLiK^] 
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Fig. 6. Spectral analysis of the Planck 70, 100, and 143 GHz maps. Columns show results computed using the three sky masks 
in Fig. 1, with, from left to right, /sky = 69.4%, 59.6%, and 39.7 %. Top row: CMB anisotropy spectra binned over a range of 
multipoles M = 40, for I > 30, with (2i + l)-weighting applied within the bin. Error bars are computed as a measure of the rms- 
power within each bin, and hence comprise both the measurement inaccuracy and cosmic variance. The grey curve is the best-fit 
Planck 6-parameter ACDM model from Planck Collaboration XVI (2014). Noise spectra computed from the half-ring-difference 
maps are shown: for the 70 GHz channel, the S/N ^ I at ^ ^ 650. Middle row: Residuals of the same power spectra with respect to 
the Planck best-fit model. Bottom row: Power ratios for the 70 vs. 100 GHz and 143 vs. 100 GHz channels of Planck. The ratios are 
calculated i by i, then binned. The error bars show the standard error of the mean for the bin. The effect of diffuse foregrounds is 
clearly seen in the changes in the 143/100 ratio with sky fraction dXi ^ 100. Bin-to-bin variations in the exact values of the ratios 
with sky fraction emphasize the importance of making such comparisons precisely. 
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Table 2. Summary of ratios of Planck 70, 100, and 143 GHz power spectra appearing in this paper. 






Spectrum Ratios 

Location 

Features 

/sky 

r Range 

70/100 

143/100 

Sect. 3.4, Fig. 6, bottom centre . . . 

No corrections 

59.6 % 

10<i< 390 

70 < ^ < 830 

1.0089 

1.0140 

1.0039 

1.0020 

Sect. 3.4, Fig. 9 . 

DSR^ correction 

59.6 % 

70 < ^ < 390 

70 < ^ < 830 

1.0080 

1.0094 

1.0045 

1.0043 

Sect. 3.4, SMICA. 

Paper VI, figure 35 


50 < ^ < 300 

300 < ^ < 700 

1.0060 

1.0075 

1.0020 

1.0020 

Sect. 4.3, Fig. 12. 

NS^ correction 

59.6 % 

70 < ^ < 390 

70 < ^ < 830 

1.0052 

1.0077 

1.0040 

1.0020 

Sect. 4.3, Fig. 13. 

DSR^-i-NS^ corrections 

59.6 % 

70 < ^ < 390 

70 < ^ < 830 

1.0043 

1.0032 

1.0046 

1.0043 

Sect. 5, Fig. 14. 

CamSpecLikelihood^ 




1.00058 


^ Discrete-source residual correction. 

^ The mask used in Paper VI, figure 35 was similar but not identical to the 39.7 % mask of Fig. 1. The differences do not affect the comparison. 
^ Near sidelobe correction, 100 and 143 GHz. 

Planck Collaboration XVI (2014). 


4.2. Nominal beam approximation 

In the 2013 analysis, both LFI and HFI performed a “nominal 
beam” calibration, i.e., we assumed that the detector response to 
the dipole can be approximated by the response of the nominal 
beam alone, which in turn is modelled as a pencil beam (for 
details see Appendix B). Clearly, if 100% of the power were 
contained in the nominal beam, the window function would fully 
account for beam effects in the reconstructed map and power 
spectrum. In reality, however, a fraction of the beam power is 
missing from the nominal beam and appears in the near and far 
sidelobes, affecting the map and power spectrum reconstruction 
in ways that depend on the level of coupling of the sidelobes 
with the dipole. Accordingly, a correction factor is applied that 
has the form (see Eq. B.12) 

Tsky ~ Tsky “ 0sky + ? (2) 

where Tsky is the true sky temperature, Tsky is the sky tem¬ 
perature estimated by the “nominal beam” calibration, 0 d = 
(Fside * ^)/(7^nominai * is the Coupling of the (near and far) 
sidelobes with the dipole, and 0 sky = (T^side * Tsky)/Tsky is a small 
term (of order 0.05 %, see Appendix C) representing the sidelobe 
coupling with all-sky sources other than the dipole (mainly CMB 
anisotropies and Galactic emission). The term 0 d is potentially 
important, since dipole signals contributing to the near sidelobes 
may bias the dipole calibration. Our current understanding of the 
value and uncertainty of the scale factors 77 = - (p^ky + 0 d) for 

LFI and HFI is discussed in detail in Appendix C. 

4.3. Key findings 

There are two key findings or conclusions from the analyses in 
Appendices B and C. 

- For LFI, a complete accounting of the corrections using the 
current full An beam model would lead to an adjustment 
of about 0.1 % in the amplitude of the released maps (i.e., 
0.2 % of the power spectra). At present this is an estimate, 
and rather than adjusting the maps we include this in our un¬ 
certainty. 


- For HFI, recent work on a hybrid beam profile, including 
data from planet measurements and GRASP^ modelling, has 
led to improvements in the beam window function correc¬ 
tion rising from 0 to 0.8, 0.8, 0.5, and 1.2 % over the range 
^ = 1 to ^ = 600, at 100, 143, 217, and 353 GHz, respec¬ 
tively. Uncertainties in these corrections have not been fully 
characterized, but are dominated by the intercalibration of 
Mars and Jupiter data and are comparable to the corrections 
themselves (see Fig. C.3) 

Figure 11 shows the corrections to the beam window func¬ 
tions at 100, 143, and 217 GHz. Figure 12 shows the effect of 
those corrections on the 70/100 and 143/100 power spectrum 
ratios, uncorrected for unresolved source residuals. There is al¬ 
most no effect on the 100/143 GHz ratio, as the dijferential beam 
window function correction between these two frequencies is 
small. The 70/100 ratio, however, is significantly closer to unity. 
In Sect. 5, we show that such a correction does not materially 
affect the 2013 cosmology results. 

Figure 13 shows the power spectrum ratios corrected for 
both the beam window functions and unresolved source resid¬ 
uals (Sect. 3.3). The average ratios over the range lb <i < 390 
are 1.0043 and 1.0046 for 70/100 and 143/100, respectively. For 
the range lb < I < 830, they are 1.0032 and 1.0043. 

For the 2014 release, we expect internal consistency and un¬ 
certainties to further improve as more detailed models of the 
beam and correction factors are included in the analysis. 

We have concentrated in this section on beam effects; how¬ 
ever, the transfer function depends also on the residuals of the 
time transfer function, measured on planets and glitches, and 
deconvolved in the time-ordered data prior to mapmaking and 
calibration. For the HFI channels, the transfer function used for 
the 2013 cosmological analysis assumes that all remaining ef¬ 
fects are contained within a 40' x 40' map of a compact scanning 
beam and corresponding effective beam. Any residuals from un¬ 
corrected time constants longer than 1 s are left in the maps, and 
will affect the dipoles and thus the absolute calibration. This has 
been investigated since the 2013 data release; time constants in 


^ Developed by TICRA (Copenhagen, DK) for analysing general re¬ 
flector antennas (http: //www. ticra .it). 
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Fig. 7. Same as the top two panels of the middle column of 
Fig. 6, but without inclusion of signal cosmic variance in the 
uncertainties. Both signal x noise and noise x noise terms are 
included. 


Fig. 8. Estimates of the residual thermal SZ and unresolved radio 
and infrared source residuals that must be removed. 

5. Likelihood analysis 

In the previous section we showed how work since release of 
the 2013 Planck results has led to an improved understanding 



the 1-3 s range have been identified and shown to be the origin 

of difficulties encountered with calibration based on the orbital Fig. 9. Same as the bottom middle panel of Fig. 6, but corrected 
dipole. The 2014 data release will include a correction of these for differences in unresolved-source residuals (see text). We have 
effects, and the absolute calibration will be carried out on the or- not tried to account for uncertainties in the foreground correction 
bital dipole. A reduction in calibration uncertainties by a factor itself; however, since the correction is small, the effect on the 
of a few can be anticipated. uncertainties would be small. 
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Deviation from the beam axis [degrees] 


Fig. 10. Radial slice through a 70 GHz beam from the GRASP 
model, illustrating the nominal beam, near sidelobe, and far side- 
lobe regions. The exact choice of angular cutoff for the nominal 
beam is different for different frequencies. 


of the beams and a small (and well within the stated uncertain¬ 
ties) revision to the near-sidelobe power in the HFI beams, which 
brings HFI and LFI into even closer agreement. In this section, 
we show that the revision in the HFI beams has little effect on 



Fig. 11. Effective beam window function corrections from 
Fig. C.3, which correct for the effect of near-sidelobe power 
missing in the HFI beams used in the 2013 results (Sect. C.2.1). 
Uncertainties are not shown here for clarity, but are shown in 
Fig. C.3, and would be large on the scale of this plot. The 
217 GHz correction is shown for illustration only. 



Fig. 12. Same as the bottom middle panel of Fig. 6, but cor¬ 
rected for the near-sidelobe power at 100 and 143 GHz that was 
not included in the 2013 results. Since the beam corrections for 
100 and 143 GHz are nearly identical, the ratio 143/100 hardly 
changes. The ratio 70/100, however, changes significantly, mov¬ 
ing towards unity. Uncertainties in the beam window function 
corrections are not included. 



0 200 400 600 800 


Fig. 13. Same as the bottom middle panel of Fig. 6, but corrected 
for both the near-sidelobe power at 100 and 143 GHz that was 
not included in the 2013 results and for unresolved source resid¬ 
uals (Sect. 3.3). Uncertainties in the beam window function cor¬ 
rections are not included. 
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I 

Fig. 14. Ratios of 100, 143, and 217 GHz power spectra calcu¬ 
lated from detector sets with the likelihood method, including 
subtraction of the best-fitting foreground model (see text) and 
correction for the best-fit relative calibration factors for individ¬ 
ual detectors. Solid symbols and lines show ratios for mask G22; 
open symbols and dotted or dashed lines show ratios for mask 
G35. The greater scatter in 217/143 for mask G35 is caused by 
CMB-foreground cross-correlations. 


cosmological parameters. To do this, we make use of the like¬ 
lihood and parameter estimation machinery described in Planck 
Collaboration XV (2014) and Planck Collaboration XVI (2014). 
For both analytical and historical reasons there are differences 
(e.g., masks, frequencies, multipole ranges) in the analyses in 
this section and in previous sections; however, as will be seen, 
the effects of the differences are accounted for straightforwardly, 
and do not affect the conclusions about parameters. 

The Planck 2013 cosmological parameter results given in 
Planck Collaboration XVI (2014) are determined for ^ > 50 
from 100, 143, and 217 GHz “detector set” data described in 
Planck Collaboration XV (2014, Table 1), by means of the 
CamSpec likelihood analysis described in the same paper that 
solves simultaneously for calibration, foreground, and beam pa¬ 
rameters. This approach allows power spectrum comparisons 
to sub-percent level precision, using only cross-spectra (as in 
Sect. 3) to avoid the need for accurate subtraction of noise in 
auto-spectra. 

In this section, we determine the ratios of the 100, 143, and 
217 GHz spectra using this approach, and compare the 143/100 
results to those found in Sect. 3. We show that the apparent dif¬ 
ference in the results from the two different approaches is easily 
accounted for by differences in the sky used, the difference be¬ 
tween the detector set data and full frequency channel data, and 
the use of individual detector recalibration factors in the detec¬ 
tor set/likelihood approach. Having established essentially exact 
correspondence between the methods, we use the likelihood ma¬ 
chinery to estimate the effect on cosmological parameters of the 
revision in the near-sidelobe power in the HFI beams. 

In Planck Collaboration XVI (2014), we used mask G45 
(/,ky = 0.45) for 100 X 100 GHz, and mask G35 (/sky = 0.37) 
for 143 X 143 GHz and 217 x 217 GHz to control diffuse fore¬ 



Fig.15. Effects on the 143/100 ratio of changes in the mask, 
choice of detectors, and detector recalibration. Solid lines in¬ 
dicate ratios calculated with mask G22; dashed lines indicate 
mask G35. Use of detector sets gives the cyan curves with recal¬ 
ibration turned on and the red curves with recalibration turned 
off. Use of full-frequency half-ring cross spectra, as in Sect. 3, 
gives the blue curves. The cyan curves are comparable to the 
green curves in Fig. 14, which also have intra-frequency calibra¬ 
tion and foreground corrections applied. The blue-dashed curve 
agrees extremely well with the blue curve in the bottom right 
panel of Fig. 6 , as it should (see text). 


grounds. However, here we are interested in precise tests of 
inter-frequency power spectrum consistency, so (as before) we 
need to compute spectra using exactly the same masks to can¬ 
cel the effects of cosmic variance from the primordial CMB. We 
have therefore recomputed all of the spectra using mask G22 
(/sky = 0.22) and mask G35, restricting the sky area to reduce 
the effects of Galactic dust emission at 143 and 217 GHz. The 
spectra are computed from means of detector set cross-spectra 
Planck Collaboration XV (2014). For each spectrum, we subtract 
the best-fitting foreground model from the Planck-l-WP-l-high^ 
solution for the base six-parameter ACDM model with param¬ 
eters as tabulated in Planck Collaboration XVI (2014), and cor¬ 
rect for the best-fit relative calibration factors of this solution. 
The convention adopted in the CamSpec likelihood fixes the cal¬ 
ibration of 143 X 143 to unity, hence calibration factors multi¬ 
ply the 100 X 100 and 217 x 217 power spectra to match the 
143 X 143 spectrum. The best fit values of these coefficients are 
cioo = 1.00058 and C 217 = 0.9974 for mask G22, both very close 
to unity and consistent with the calibration differences between 
individual detectors at the same frequency (see Table 3 of Planck 
Collaboration XV 2014). The results are shown in Fig. 14. 

The 143/100 ratio given by the dashed green line can be com¬ 
pared with the bottom right panel of Fig. 6 , which is based on 
40 % of the sky, nearly the same as mask G35. As expected, they 
are not identical; Fig. 15 explains the differences. In Fig. 15, 
pairs of curves in the same colour show the difference between 
mask G22 and mask G35, as labelled. The cyan curves can be 
compared to the green curves in Fig. 14, which also have inter¬ 
frequency calibration and foreground corrections applied. The 
red curves show the effect of turning off detector-by-detector 
intercalibration. The blue curves show the effect of switching 
from detector sets to full-frequency half-ring cross spectra (as 
in Sect. 3). The progression from solid cyan to dashed blue in 
Fig. 15 shows the relationship between the PLA map-based re¬ 
sults and the detector-set/likelihood results. As used in the like- 
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lihood analysis (Planck Collaboration XV 2014), the 143/100 
ratio is 1.00058 over the full i range used in the likelihood anal¬ 
ysis, compared to the ratios between 1.0039 and 1.0046 seen in 
Table 2 over 70 < ^ < 390 for Figs. 6, 9, 12, or 13. However, us¬ 
ing mask G35 (/sky = 0.37), using half-ring cross-spectra of full- 
frequency detector sets, and turning off unresolved-source resid¬ 
ual and detector-by-detector intercalibration factors, changes the 
ratio over 6Q < i < 390 to 1.0033, in good agreement with the 
1.0039 calculated for the 143/100 comparison in the bottom right 
panel of Fig. 6. 


This agreement extends to the detailed shapes of the two 
curves (blue in the bottom right panel of Fig. 6 and blue-dashed 
in Fig. 15) as well. This is necessarily the case, since they are 
both cross-spectra of half-ring frequency maps, without cor¬ 
rections for unresolved-source residuals, and using the “2013” 
beams. The only difference in the data comes from the masks 
used, which are the GAL040 mask (/sky = 39.7 %) and mask 
G35 (/sky = 37 %), respectively. This agreement is nevertheless 
reassuring in showing that the differences in spectral ratios be¬ 
tween the PLA map-based approach and the detector set likeli¬ 
hood approach are well-understood, and disappear for common 
data and masks. 

We can now turn to the question of whether the small revi¬ 
sion to the HFI beams affects cosmological parameters. A full re¬ 
vised beam analysis at the detector level that includes the 0.1 % 
power in near sidelobes not taken into account directly in the 
100 and 143 GHz beams in 2013 (Sect. 4) has not yet been 
completed; however, for an indicative test, we rescaled the aver¬ 
aged cross-spectra appearing in the likelihood by functions cor¬ 
responding to the new beam shapes for the ^-ranges for which 
they have been calculated (presently up to ^ = 2000). Where 
necessary the shapes were extrapolated as being flat up to higher 
i. The 143 X 217 spectrum was rescaled by the geometric mean 
of the 143 and 217 rescalings. Then we performed a Monte 
Carlo Markov Chain (MCMC) analysis for the base ACDM 
model for the modified “high-7” ’ likelihood with an unmodified 
low-7 Planck likelihood and WMAP low-7 polarized likelihood 
(“WP”). To see any change in the beam error behaviour, we 
choose to sample explicitly over all twenty of the eigenmode am¬ 
plitudes, rather than sampling over one and marginalizing over 
the other nineteen, as we did in the parameters paper (Planck 
Collaboration XVI 2014). 

The results are indicated in Figs. 16 and 17, showing a selec¬ 
tion of cosmological parameters and the beam eigenmode am¬ 
plitudes, respectively. As expected, we see a boost in the power 
spectrum amplitude, resulting in a change to the cosmological 
amplitude at about the 1 cr level. However, the largest shift in any 
other cosmological parameter is 0.3 cr. The uncertainty in the 
beam window function is described by a small number of eigen- 
modes in multipole space and their covariance matrix (Planck 
Collaboration VII2014). The posteriors for the first beam eigen- 
modes for the 100, 143, and 217 effective spectra shift notice¬ 
ably; others are practically unchanged. The beams used here are 
preliminary and the beam eigenmodes have not been generated 
self-consistently to match the beam calibration pipeline. No ad¬ 
justment was made in the calibration of the low-7 likelihood. 
Nevertheless, from the results presented here, we can anticipate 
that the 2014 revisions to the beams will affect the overall cal¬ 
ibration of the spectra, but will have little other impact on cos¬ 
mology. 
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Fig. 16. Changes in cosmological parameters from the inclusion 
of the near sidelobe power discussed in the text. The black curves 
are the 2013 results for Planck plus the low-7 WMAP polariza¬ 
tion (WP). The red curves are for Planck+WV using the revised 
HFI beams. The shifts in the posteriors are all less than 0.3(T ex¬ 
cept for the cosmological amplitude As and parameters related 
to it, as expected. 


6. Comparison of Planck and WMAP 

Planck and WMAP have both produced sky maps with excellent 
large-scale stability, as demonstrated by many null tests both 
internal to the data and external. In this section, we compare 
Planck and WMAP measurements in several different ways. In 
Sect. 6.1, we compare power spectra calculated from 70 and 
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Fig. 17. Changes in the beam eigenmode coefficients from the inclusion of the near-sidelobe power now established but not included 
in the processing for the 2013 results. Black curves are for Planck+^V\ red curves are for Planck+^V with the revised beams. The 
superscript indicates the effective spectum (one to four for 100, 143, 217, and 143 x 217 respectively) while the subscript indicates 
eigenmode number. 


100 GHz Planck maps available in the PL A, and from V- and 
W-band yearly maps in the WMAP9 data release. In Sec. 6.2, 
we perform a likelihood analysis similar to that in Sect. 5 and 
in Planck Collaboration XVI (2014), and show that the differ¬ 
ences between the map-based and likelihood analyses are well- 
understood. In Sect. 6.3 we assess the results in the context of 
the uncertainties for the two experiments. 


6.1. Map and power spectrum analysis 

The WMAP9 data release includes Aside = 1024 yearly sky maps 
from individual differential assemblies (DAs), both corrected 
for foregrounds and uncorrected, as well as Aside = 512 fre¬ 
quency maps. WMAP uses somewhat different sky masks than 
Planck. In Sect. 3 we emphasized the importance of using ex¬ 
actly the same masks in comparing results. Accordingly, for 
PlanckjWMAP comparisons we construct a joint mask, taking 
the union of the Planck GAL060 mask used in Sects. 3 and 
4, the WMAP KQ85 mask, which imposes larger cuts for ra¬ 
dio sources and some galaxy clusters, as required by the poorer 
angular resolution of WMAP, and the Planck joint 143, 100, 
70 GHz point source mask. Fig. 18 shows the mask, which leaves 
/sky = 56.7 % of the sky available for spectral analysis.. 


We use the same spectrum estimation procedure as in Sect. 3, 
evaluating the relevant cross-spectra, correcting for the mask 
with the appropriate kernel, and dividing out the relevant beam 
response and pixel smoothing functions. As the mask is different 
from the one used for Planck-only comparisons, so is the mask- 
correction kernel. All maps are analysed using the same mask. 

For WMAP, there are nine yearly sky maps for each differ¬ 
ential assembly VI, V2, Wl, W2, W3, and W4, at Aside = 1024. 
Because the WMAP V band and the Planck 70 GHz band are so 
close in frequency, as are W band and 100 GHz, we use maps 
not corrected for foregrounds for the comparison. All possible 
cross spectra from the yearly maps and differential assemblies 
are computed (630 at W band, 153 at V), and corrected for the 
mask, beam (using WMAP beam response functions, different 
for each differential assembly), and pixel-smoothing. The cor¬ 
rected spectra are averaged, and the error on the mean is com¬ 
puted for each Q. These average differential-assembly spectra 
are then co-added with inverse noise weighting to form one 
V band and one W band spectrum. These are binned (^min = 30, 
Al = 40), and rms errors in the bin values are computed. The 
resulting spectra are shown in Fig. 19. 

The 70, 100, and 143 GHz Planck spectra and spectral ratios 
in Figs. 19-21 are determined as before, but using the new mask. 
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Fig. 18. Planck /sky ~ 60 % Galactic mask in yellow and WMAP 
KQ75 at A^side = 1024 in red. The mask used for compara¬ 
tive spectral analysis of the Planck 70, 100, and 143 GHz, and 
WMAP nine-year V- and W-band sky maps is the union of 
the two. The joint Planck 70, 100, and 143 GHz point source 
mask is also used, exactly as before. The final sky fraction is 
/sky = 56.7 %. The Planck mask is degraded to the pixel resolu¬ 
tion A^side = 1024, at which the yearly WMAP individual differ¬ 
ential assembly maps are available. 


starting from the 70 GHz A^side = 1024 half-ring PLA maps and 
the 100 and 143 GHz Aside = 2048 maps degraded to Aside = 
1024. Thus all spectra are evaluated with the identical mask, at 
the same resolution. Spectral binning and the estimation of rms 
bin errors proceed in exactly the same way as for the WMAP 
spectra and for previous Planck-only comparisons. 

Figure 19 compares the Planck 70 GHz power spectrum with 
the WMAP V-band spectrum, and the Planck 100 GHz power 
spectrum with the WMAP W-band spectrum. The Planck 2013 
best-fit model is shown for comparison. The Planck and WMAP9 
spectra disagree noticeably in the ^-range of the first two peaks. 
Ratios of spectra in Fig. 20 show this disagreement directly. In 
Figs. 20-22 the 70/100 and 143/100 ratios are the same as in 
Sects. 3 and 4, except for the small change in the mask. 

Figure 21 is the same as Fig. 20, but with the Planck 70/100 
and 143/100 ratios corrected for the missing near sidelobe power 
in the 100 and 143 GHz channels, discussed in Sect. 4. 

Figure 22 is the same as Fig. 21, but with all spectra addi¬ 
tionally corrected for residual unresolved sources, as described 
in Sect. 3.3. Mean values of the ratios over specified multipole 
ranges are given in Table 3. 

6.2. Likelihood analysis 

The likelihood analysis here is slightly different from the anal¬ 
ysis presented in Sect. 5. First, we created a point source mask 
by concatenating the WMAP/70/100/143/217 point source cata¬ 
logues. We present here only the results comparing WMAP, LFI 
70 GHz, and HFI100 GHz. Restricting the frequencies to a range 
close to the diffuse foreground minimum has the advantage that 
we can increase the sky area used. We therefore present results 
for mask G56 (unapodized), which leaves 56 % of the sky. This 
mask, combined with the point source mask, is degraded in res¬ 
olution from Aside = 2048 to the natural Aside for WMAP (512) 
and LFI (1024). 

Beam-corrected spectra for the WMAP V, W, and V-hW 
bands, and for the LFI 70 GHz bands were computed. Errors 
on these spectra were estimated from numerical simulations. At 
70 GHz we have three maps of independent subsets of detec¬ 
tors, therefore we can estimate three pseudo-spectra by cross- 
correlating them by pairs. Pseudo-spectra are then mask- and 


beam-deconvolved. The final 70 GHz spectrum is obtained as a 
noise-weighted average of the cross-spectra. Noise variances are 
estimated from anisotropic, coloured noise MC maps for each 
set of detectors. We use the FEBeCoP effective beam window 
functions for each subset of detectors (the three pairs 18-23, 19- 
22, and 20-21 as defined in Planck Collaboration II 2014). In 
the present plots we do not show beam uncertainties, which are 
bounded to be ^ 0.2 % in the multipole range consid¬ 

ered here (Planck Collaboration VI 2014). 

We also derive power spectra for the WMAP V and W bands, 
using coadded nine year maps per DA, for which no foreground 
cleaning has been attempted. There are two DA maps for V-band 
and four for W-band. The V-band spectrum is obtained by cross- 
correlating the two available maps, whereas the W spectrum is 
the noise-weighted average of the six spectra derived by cor¬ 
relating pairs of the four DA maps. We have produced Monte 
Carlo simulations of noise in order to assess the error bars of 
the WMAP spectra. We generate noise maps according to the 
pixel noise values provided by the WMAP team rescaled by the 
number of observations per pixel. Beam transfer functions per 
DA are those provided by the WMAP team. In the present error 
budget we did not include beam uncertainties, which would be 
AB^IB^ ^ 0.4 % for V and 0.5 % for W, over £ < 400. 

For LFI and WMAP, we subtract an unresolved thermal SZ 
template with the amplitude derived from the CamSpec like¬ 
lihood analysis, which we extrapolate to the central LFI and 
WMAP frequencies. We also subtract a Poisson point source 
term with an amplitude chosen to minimize the variance be¬ 
tween the measured spectra and the best-fit theoretical model at 
high multipoles. We then determine calibration factors cx for the 
WMAP and 70 GHz LFI spectrum relative to the 100 GHz HFI 
spectrum by minimizing 


A 




(3) 


until convergence, at each iteration determining cr\, an “excess 
scatter” term, by requiring that the reduced be unity. (The 
excess scatter comes primarily from foreground-CMB correla¬ 
tions, as discussed in Planck Collaboration XV 2014.) The sum 
extends over bins with central multipole in the range 50 < ^ < 
400. The index X denotes the spectrum (70, V, W, or V-i-W), and 
CTxb is the noise contribution to the error in band b. The spec¬ 
tra D in Eq. (3) are foreground-corrected. Calibration factors 
Cx and resulting spectra for 70 GHz, V, W, and V-rW-band rel¬ 
ative to 100 GHz are given in Fig. 23 and Table 3. The value of 
C 70 = 0.994 seen in Fig. 23 is entirely consistent with the 70/100 
ratios given in Sect. 5, taking into account mask and dataset vari¬ 
ations as discussed in Sect. 5. A calibration factor for V relative 
to 70 GHz calculated the same way (D^^^ is replaced by in 
Eq. 3) is also included in Table 3. 

To compare the likelihood results with the map-based ones, 
we need to identify the closest cases for which results are given. 
No correction for near sidelobes at 100 GHz (Sect. 4.3) was in¬ 
cluded in the likelihood results, but corrections for residual unre¬ 
solved sources were. On the other hand, from the results shown 
in Table 3, residual-source corrections make negligible differ¬ 
ence. In the likelihood approach, £ values beyond about 300 will 
be significantly down-weighted due to increased variance. Thus 
the most reasonable comparison is between the ratios from the 
map-based approach in the £ range 110 < £ < 310 with no cor¬ 
rections, and the corresponding calibration factors from the like¬ 
lihood approach. Values in Table 3 show good agreement be¬ 
tween the two analyses. We have furthermore noticed that using 
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Fig. 19. Left: WMAP V band compared to Planck 70 GHz and best-fit model. Right: WMAP W band compared to Planck 100 GHz 
and best-fit model. The joint PlanckjYJ^S sky mask -i- Planck point source mask (/sky = 56.7 %; see Fig. 18) is used. Because the 
frequencies are so close, no corrections for foregrounds are made. 

Table 3. Summary of ratios of Planck 70 and 100 GHz and WMAP V- and W-band power spectra appearing in this paper. 







Spectrum Ratios 

Location 

Features 

/sky 

i Range 

70/V 

100/W 100/(V-fW) 

Sect. 6.1, Fig. 20. 

No corrections 

56.7 % 

70 < ^ < 390 
110<^<310 

0.983 

0.981 

0.979 

0.977 

Sect. 6.1, Fig. 21. 

Near-sidelobe (NS) correction, 100 GHz 

56.7 % 

70 < ^ < 390 
110<^<310 

0.983 

0.981 

0.983 

0.981 

Sect. 6.1, Fig. 22. 

Discrete-source residual-i-NS corrections 

56.7 % 

70 < ^ < 390 
110<^<310 

0.983 

0.981 

0.983 

0.981 

Sect. 6.2. 

Likelihood “calibration factor”; WMAP full-mission 
maps at A^side = 512 rather than yearly maps at 
A^side = 1024; Planck detector sets (see Sect. 5) 
rather than full-frequency half-ring data. 

56.7 % 

50 < ^ < 400 

0.978 

0.976 0.974 


the WMAP full-mission maps at resolution A^side = 512 rather 
than the yearly maps at A^side = 1024, and Planck detector sets 
rather than full-frequency half-ring maps, can lower the spectral 
ratio by 0.2 % in the multipole range 110 < ^ < 310. 

6.3. Assessment 

Both the direct comparison of power spectra from the 2013 re¬ 
sults maps and the likelihood analysis show a discrepancy be¬ 
tween Planck and WMAP across the region of the first peak, 
where the S/N for WMAP is good. This difference is about 2 % in 
power, corresponding to 1 % in the maps. These numbers quan¬ 
tify what can be seen by eye in Fig. 19. The result is roughly 
the same for comparisons of 70 GHz with V-band and 100 GHz 
with W-band, where the frequency differences are small enough 
to rule out foregrounds as the cause. There is some variation of 
the ratio with i (perhaps seen more easily in Fig. 23), suggest¬ 


ing that the cause is not simply a result of calibration errors, but 
neither does the shape correspond obviously to what would be 
expected from missing power in beams. 

The 70/V and 100/W ratios differ from unity by more than 
expected from the uncertainties in absolute calibration deter¬ 
mined for Planck and WMAP. Calibration of WMAP9 is based 
on the orbital dipole (i.e., the modulation of the solar dipole due 
to the Earth’s orbital motion around the Sun), with overall cali¬ 
bration uncertainty 0.2 % (Bennett et al. 2013). The absolute cal¬ 
ibration of the Planck 2013 results is based on the solar dipole 
(Sect. 1), assuming the WMAP5 value of (369.0 + 0.9)kms“\ 
where the 0.24 % uncertainty includes the 0.2 % absolute cal¬ 
ibration uncertainty (Hinshaw et al. 2009). The overall cali¬ 
bration uncertainty is 0.62 % in the 70 GHz maps and 0.54 % 
in the 100 and 143 GHz maps (Planck Collaboration I 2014, 
Table 6). When comparing Planck and WMAP calibrated maps, 
one must remove from these uncertainties the 0.2 % contribu- 
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Fig. 20. Ratios of power spectra for Planck and WMAP over the joint Planck/KQ75-\-Planck point source mask with /sky = 56.7 %. 
The Planck 70/100 and 143/100 ratios can be compared to the bottom middle panel in Fig. 6 for /sky = 59.6%. Here we limit 
the horizontal scale because the WMAP noise beyond I = 400 makes the ratios uninformative. The general shape of the Planck 
ratios is the same; however, it is clear (as it was in Fig. 6) that changes in the sky fraction change the ratios significantly within the 
uncertainties. This underscores the importance of strict equality of all factors in such comparisons. 


tion of the WMAP absolute calibration uncertainty to the WMAP 
dipole, which thus affects both experiments equally. The remain¬ 
ing 0.14 % uncertainty in the WMAP dipole (mostly due to fore¬ 
grounds; Hinshaw et al. 2009) affects Planck but not WMAP, 
because its absolute calibration is from the orbital dipole. In 
the planned 2014 release, the Planck absolute calibration will 
be based on the orbital dipole, bypassing uncertainties in the 
solar dipole. At the power spectrum level for comparison with 
WMAP, then, the Planck uncertainty would be between 2 x 
(0.54 - 0.2) % = 0.68 % and 2 x (0.62 - 0.2) % = 0.84 %, and the 
WMAP uncertainty 0.4 %. The power spectrum ratios in Table 3 
from Fig. 21 and Sect. 6.2 then represent a 1.5-2(T difference. 

As the primary calibration reference used by Planck in the 
2013 results is the WMAP solar dipole, the inconsistency be¬ 
tween Planck and WMAP is unlikely to be the result of sim¬ 
ple calibration factors. Reinforcing this conclusion is the fact 
that the intercalibration comparison given in Fig. 35 of Planck 
Collaboration VI (2014) for CMB anisotropies shows agree¬ 
ment between channels to better than 0.5 % over the range 70- 
217 GHz, and 1 % over all channels from 44 to 353 GHz, us¬ 
ing 143 GHz as a reference. Problems with transfer functions are 
more likely to be the cause. The larger deviations at higher mul¬ 
tipoles in the Planck intercalibration comparison just referred to 
also point towards transfer function problems. 

Comparisons between WMAP- and LFI-derived bright¬ 
ness temperatures of Jupiter presented in figure 21 of Planck 
Collaboration V (2014) provide a potentially useful clue: the two 
instruments seem to agree (for a simple linear spectral model for 
Jupiter) within 1 % , and much better than this at 70 GHz. The 
LFI points are derived from just two Jupiter transits. We now 
have eight transits in hand. If the full analysis of all the Jupiter 


crossings shows similar consistency with WMAP brightness tem¬ 
peratures, subtle effects from beam asymmetries and near side- 
lobes in the two instruments, rather than gross errors in main 
beam solid angles, may be the cause. 

At present, we do not have an explanation for the 2% cal¬ 
ibration difference between WMAP and Planck. The differences 
between WMAP and Planck are primarily multiplicative in the 
power spectra, and so have little impact on cosmological param¬ 
eters other than on the amplitude of the primordial spectrum A^. 
Figure 23 shows that after the spectra are rescaled by multiplica¬ 
tive factors, there is excellent point-to-point agreement between 
the LFI, HFI, and WMAP power spectra. Appendix A of Planck 
Collaboration XVI (2014) also shows explicitly that by restrict¬ 
ing the Planck 2013 likelihood to ^max = 1000, we recover al¬ 
most identical cosmological parameters to those from WMAP 
apart from a small shift in As reflecting the calibration differ¬ 
ences between the two experiments. Any shape differences be¬ 
tween the power spectra of the two experiments therefore have 
little impact on cosmological parameters. 


7. Conclusions 

Data consistency in CMB experiments is an important indica¬ 
tor of the validity of the results derived from these data. We have 
compared calibrated data from different Planck channels in three 
ways: (1) comparison of power spectra calculated on masked, 
publicly-released maps at 70, 100, and 143 GHz, including cor¬ 
rections for residual unresolved sources; (2) comparison of cal¬ 
ibration factors determined to minimize differences between the 
CMB in the maps; and (3) comparison of power spectra and cal¬ 
ibration factors obtained in the likelihood analysis used to de- 
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Fig. 21. Same as Fig. 20, but the Planck 70/100 and 143/100 ratios are corrected for beam power at 100 and 143 GHz that was not 
included in the effective beam window function used in the 2013 results. 



Fig. 22. Same as Fig. 20, but including corrections for both Planck beams and for Planck and WMAP unresolved-source residuals. 


termine cosmological parameters. We have emphasized particu¬ 
larly the comparison of 70 GHz and 100 GHz, as these frequen¬ 
cies come from different instruments and are analysed indepen¬ 
dently, and are near the diffuse foreground minimum. The inter¬ 
nal consistency of the Planck data is remarkable (relative cal¬ 
ibration inconsistencies between CMB channels of only a few 
tenths of a percent), and is at a level consistent with the uncer¬ 


tainties based on the 2013 solar dipole calibration scheme. In 
the future, calibration based on the orbital dipole should lead to 
more accurate absolute calibration. 

The estimate of the contribution of the near sidelobes to the 
HFI beam solid angles in March 2013 implicitly assumed az¬ 
imuthal symmetry of the near sidelobes. Subsequent analyses of 
the beam and calibration procedures suggest that this assump- 
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Fig. 23. Foreground-corrected spectra computed using mask 
G56, with calibration factors Cx from Eq. 3 applied to 70 GHz, 
V, W, and V-i-W. The solid red lines in the upper panels show the 
best-fit model CMB spectrum (Planck Collaboration XVI2014). 
The lower panels show the residuals with respect to this model. 
The error bars on the LFI and WMAP points show noise errors 
only. 


tion underestimated near-sidelobe power by as much as 0.1 %. 
Corrections for this missing power result in improved agreement 
between 70 and 100 GHz, measured by ratios of power spectra, 
at a level of typically 0.3 %, well below the levels of uncertainty 
estimated for LFI and HFI calibration, and therefore having little 
effect on the cosmological parameters determined from the 2013 
data. 

Similar analysis applied to WMAP data shows a roughly 2 % 
difference between Planck and WMAP over the region of the first 
acoustic peak. At present, the explanation for this discrepancy 
with WMAP is not known. 

Precision calibration is an ongoing challenge. Future analy¬ 
ses will incorporate the lessons of this study of consistency, fur¬ 
ther reducing the size of uncertainties and errors on the scientific 
results from Planck. 
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Appendix A: Difference maps 

Figure A. 1 shows the frequency maps and difference maps used 
to calculate the rms values given in Table 1 . Frequency maps are 
half-ring sums. Difference maps are half-ring differences. The 
maps have been processed as described in Sect. 2 to retain struc¬ 
ture between angular scales of 15' and 8°. Figure 4 gives the 
corresponding histograms. 
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Fig. A.l. Frequency maps and difference maps using the GAL060 mask. Rms values are given in Table 1; histograms are shown in 
Fig. 4. In the top 3x3 group, the frequency maps are on the diagonal. The three frequency difference maps are in the upper right 
block. The noise difference maps are in the lower left block. The bottom row gives the single-frequency noise maps. The caption of 
Fig. 4 describes characteristics of the histograms that can also be seen quite clearly in these maps. The three upper right maps are 
the same as the three in Fig. 3, but with structure on scales larger than 8° removed, as described in the text. 
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Appendix B: Nominal beam calibration: definitions, 
approximations, and rescaling factors 


a fraction (1 - /si) of the antenna gain is contained in the nomi¬ 
nal beam, so we have 


In the Planck scanning strategy, for each pointing period, the 
spin axis remains fixed and the output voltage sees a spin- 
synchronous modulation as the telescope points at direction x 
along a scan circle. The output signal (in arbitrary units) mea¬ 
sured by a given Planck detector, either LFI or HFI, can be writ¬ 
ten as 

Vout(x) = g[p * (Tsky + £>)](x) + M{p), (B.l) 

where P is the angular response (“beam”), normalized to unit 
solid angle (£^PdQ = 1), D is the dipole signal (including 
the solar and orbital terms), Tsky is the sky brightness temper¬ 
ature (CMB and foregrounds, i.e., everything other than the 
dipole), and M is a small monopole term assumed to be con¬ 
stant for each pointing period p. We use the convolution notation 
A*P = £^APdQ. 

For each pointing period, the output voltage measures a spin- 
synchronous modulation dominated by the dipole, which is used 
to recover a gain calibration factor. The dipole term, F * Z), is 
dominant (by two orders of magnitude) once strong foreground 
regions are masked, so for calibration purposes we have Vout = 
G(P * D). For each pointing period, calibration is performed by 
fitting Vout to a model of the beam-convolved dipole: 


(P * Z))uiodel 


^nominal * ~ (1 - /sl)Fpencil ^ D = (1 - fsl)D, (B.6) 

where /si is the solid angle fraction of the sidelobes. For the 
70, 100, and 143 GHz channels, /si 0.7, 0.5, and 0.5% (see 
Sect. 4.1). In this formalism, we assume that the time response 
of the bolometers has been accurately deconvolved; thus no ad¬ 
ditional time constant corrections are needed. This approxima¬ 
tion is true to high accuracy (< 0.2 %), since the dipole changes 
very little within the nominal beam solid angle. Rearranging the 
terms, we find 

G = G(1-/si)(1+0d), (B.7) 

where /d = (^side * ^)/(^nominal * ^) represents the coupling of 
the dipole with sidelobes and depends on the detailed structure 
of the sidelobes and the scanning strategy. Equation B.7 gives 
the correction factor that one should apply to recover the proper 
gain factor if one models the beam as a pencil beam. 

B.2. Rescaling factor 

To obtain the reconstructed sky signal after calibration, rsky(x), 
the dipole term needs to be subtracted from the timelines 
(Eq. B.l). Since we are interested in the signal from the sky 
entering the nominal beam, the coupling of Tsky with sidelobes 
represents a spurious additive signal. So we write the timeline as 


where G is the estimated gain, while G = Vout I (P * D) is the true 
gain. For the 2013 release, we rely on the a priori knowledge 
of the solar dipole from WMAP (Hinshaw et al. 2009). Relative 
changes of the gain G(p) between sky circles contribute to the 
monopole term M(p), and are removed by destriping. Assuming 
an error-free knowledge of the true dipole and of the full beam 
pattern, then Vout = (P ^ T))niodei and G = G, so one would 
recover the true value of the gain, preserving the entire power 
entering the beam. However, full beam convolution was not im¬ 
plemented in the LFI and HFI pipelines for the 2013 release, and 
we adopted a “nominal beam approximation” approach. 


B.1. Nominal beam approximation 


We approximate the beam pattern with a “pencil beam,” i.e., a 
normalized delta-function: 


(P * Z))uiodel 




^pencil^model^f^ ~ ^pencil * T)uiodel* (B.3) 


Neglecting intrinsic errors in the model dipole (Dmodei ~ ^), the 
estimated gain is 


G = G 


(F*Z)) 

(P pencil * P^) 


(B.4) 


We can write P as the sum of two terms, P = ^nominal + ^side^ 
where ^nominal represents the contribution of the nominal beam 
(defined as the portion of the beam defining the window func¬ 
tion, see Sect. 4) and Fside is the contribution from the (near and 
far) sidelobes. We have 


^ ^ (^nominal * ^) (Pside * ^) 

(^pencil * P^) 


(B.5) 


The convolution of the nominal beam to the dipole is nearly 
identical to an ideal pencil beam convolution, except that only 


Vout(-^) = G(Fnominal * Psky + Pside * 7^sky)(-^) + M(p). (B.8) 


In the ideal case of an error-free, full-convolution model, G = G, 
one gets 



P nominal * P sky P side * P sky • 


(B.9) 


For angular scales larger than the nominal beam, a relation sim¬ 
ilar to Eq. (B.6) holds for Tsky, i.e., Tnominal * Tsky ~ (1 — /sl)Tsky. 
With this and Eq. (B.9), we obtain 


sky 


= T. 


sky 


1 /sky 
l-/sl ’ 


(B.IO) 


where /sky = (Tside * Tsky)/Tsky represents the coupling of the 
sidelobes with the sky signal (other than the dipole), i.e., the 
straylight pickup. 

As mentioned, however, in the 2013 release both LFI and 
HFI used the nominal beam approximation. In this case, the es¬ 
timated gain is given by Eq. (B.7), and we have 


Lout (Tnominal * Tsky) + (Tside * Tsky) 

G (1 - /si)(l +/d) 


(B.ll) 


Solving for the true sky temperature we find 


' sky 


- sky 


/sky 

(l-/sl) 


+ /d 


Tsky (l - /sky + /d) • (B.12) 


This expression does not contain the term /si to first order; how¬ 
ever, the (smaller) correction terms /d and /sky need to be evalu¬ 
ated to derive the proper normalization of the temperature maps 
and of the power spectra. 

In summary, to evaluate the true sky temperature from the 
estimated one, some correction factors need to be evaluated to 
properly account for the beam. These correction terms depend on 
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the model assumptions. For a full convolution model, the rescal¬ 
ing factor is 77 = (1 - 0 sky)(l - /si)- For a pencil beam approxima¬ 
tion, currently applied by both Planck instruments, the rescaling 
factor is 77 = (1 - 0sky + /d)- 

If one were to apply a correction, a rescaling factor r] could 
be directly applied to the temperature map to give a properly 
calibrated map. Alternatively, the window function should be 
rescaled by 77 “^. 

B.3. Semi-analytical expression for0D 

A useful expression for 0 d can be obtained as follows. The 
dipole seen through any particular line of sight (LOS) will be 
proportional to the dot product of the dipole vector with that of 
the line of sight: 

D = Do [sin 0LOS cos (Olos - ^d) + cos 0 d sin 0los] , (B.13) 

where Dq = 3.35 mK, the amplitude of the dipole, 0 d and Od 
are the polar and azimuthal angles of the direction of the dipole 
maximum on the sky^, and 0 los and Olos are the same for the 
direction of the nominal beam observation. For each ring, choos¬ 
ing a coordinate system for which the z-axis is aligned with the 
satellite spin axis, and around which both the LOS and the dipole 
axis will rotate, ensures that the last term in Eq. (B.13) is con¬ 
stant. Since a constant is removed for each stable pointing pe¬ 
riod (Planck Collaboration VIII 2014), we ignore this and write 

D = Dq sin 0LOS cos (Olos “ Od) . (B.14) 

Specializing, when the line of sight (LOS) is the nominal beam, 
the signal seen will be 

^main ~ Dq sin Omain COS (Oj^ain “ Oj)) • (B.15) 

Similarly, when we observe a point on the sky through a side- 
lobe, we have 

Ai = Do sin 0si cos (Omain + AOsi - Od) , (B. 16) 

where AOsi = Osi-Omain is the difference in azimuthal angles be¬ 
tween the nominal beam and the sidelobe of interest. This yields: 

Al = Do sin 0sl cos (AOsi) cos (Omain “ Od) 

-Do sin 0si sin (AOsi) sin (Omain - Od) . (B. 17) 

Because Planck calibrates by “locking in” on, or fitting to, 
a sinusoidal signal of the form cos (Omain - Od) (Eq. B.15), the 
second term in Eq. (B.17) is irrelevant, and the far sidelobes will 
make a contribution to the gain according to 

^side *D (0sl, ®sl) sin ©si cos (AOsi) 

^-.(B.18) 

^nominal * D Ajl Sm 0 main 

Integrations over GRASP models of the HFI far sidelobes in¬ 
dicate that 0 D ~ /fsl/3. 


^ Note that these angles should include the change to the dipole in¬ 
duced by the satellite motion with respect to the Sun, in addition to that 
induced by the Sun’s motion with respect to the CMB. We have ne¬ 
glected such boosting effects so far - see Planck Collaboration XXVII 
(2014). 


LFI18S, monochromatic beam 



Fig. C.l. 0D for a simulation using a 70 GHz beam. 


Appendix C: Assessment of LFI and HFI with 
respect to beams and caiibration 

C.1. LFI 

C.1.1. Evaluation of 0sky 

The term 0sky turns out to be small, but it is difficult to quantify 
because it requires a full convolution of a sky model with the 
full beam, and it is frequency dependent. We have performed a 
simplified simulation for a 70 GHz channel (LFI18M), and find 
values for 0sky ranging between 0.05 % and 0.2 % throughout a 
full survey. More realistic simulations have been produced re¬ 
cently on an LFI 30 GHz channel (LFI27M), for both Survey 1 
and Survey 2. We find: 

Survey 1, (/sky) = 0.05 %, 0.01 % < /sky < 0.23 %; 

Survey 2, (/sky) = 0.06 %, 0.01 % < /sky < 1.11 %. 
(Differences between odd and even Surveys are to be ex¬ 
pected, given the differences in scanning described in Planck 
Collaboration I (2014) section 4.1.) These results are consistent 
with our preliminary findings, and indicate a typical value of 
/sky ~ 0-05 % at all Planck frequencies. 

C.1.2. Evaluation of /d 

We can estimate /d by taking average values of the convolu¬ 
tions involved in its definition in Eq. (B.12). The denominator 
(^nominal ^ D) = 3.35 mK is the dipole amplitude. To estimate 
(Fside * D), we have computed the convolution of the LFI side¬ 
lobes, modelled with GRASP, with the dipole. 

Note that /d contains the effects of both near and far side¬ 
lobes; however, it is the dipole modulation through the near 
sidelobes that contributes most to its value. Destriping removes 
the contribution from sidelobes nearly aligned with the space¬ 
craft spin axis (such as the “primary spillover,” see Planck 
Collaboration I 2014), which produces a nearly constant offset 
for each pointing period. After destriping, we find that (Pside * D) 
is approximately 3 pK at 70 GHz and 5 pK at 30 GHz, leading to 
(/d) ~ 0.1 % and (/d) ~ 0.15 %, respectively. 

Figure C.l shows the pixel histogram for a simulation where 
the value of /d was computed in every sky direction for an LFI 
70 GHz channel. The median value, /d ~ 0.16%, is in good 
agreement with our first-order estimate. 
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In conclusion, the correction terms 0 d and 0sky are of the 
same order (0.15 % and 0.05 %, respectively) and tend to can¬ 
cel each other in Eq. (B.12), resulting in a net effect of 0.1 % 
on the LFI calibrated map. For the 2013 release, we have not 
carried out a systematic simulation for all channels; rather we 
have assumed a conservative residual uncertainty of 0.2% in 
the FFI calibrated maps due to sidelobe convolution effects 
(Planck Collaboration V 2014, Table 8). More detailed analysis 
is planned for the 2014 release. 

C.2. HFI 

The effective beam window function, including the effects of 
both calibration and map reconstruction, is based on measure¬ 
ments of the scanning beam (Planck Collaboration VII 2014). 
These are obtained from Mars observations after deconvolution 
of the time response of the detection chain. The effective beam 
window function is used to correct the cross-power spectra used 
in the cosmological analysis. Figure 10 of Planck Collaboration 
VII (2014) shows this scanning beam, and provides estimates 
of the near sidelobe contribution to the solid angle between 0.1 
and 0.2 % at 100 and 143 GHz, respectively. In the subsections 
below, we estimate the impact of the approximations used in es¬ 
timating the nominal beam and near sidelobe region (up to 5° 
from the beam centre) through the use of Saturn and Jupiter ob¬ 
servations. We estimate the bias introduced by these approxi¬ 
mations, as well as the bias in the calibration factor due to the 
far sidelobes. As previously stated, we make no attempt to re¬ 
calibrate the 2013 data release, but only show that the FFI-HFI 
differences at intermediate multipoles are understood within the 
present uncertainties. 

C.2.1. Nominal beam and near sidelobes 

The nominal beam for HFI is measured to roughly 15' with Mars 
observations. This provides most of the information needed for 
the high-^ analysis used to obtain the 2013 Planck cosmological 
results. Jupiter observations show shoulders extending to 20' for 
all channels (see fogire 10 of Planck Collaboration VII 2014). 
Features in this region are small and difficult to measure. Their 
integrated contribution, however, is not negligible at the level 
of a few tenths of a percent. Table 2 and figure 10 of Planck 
Collaboration VII (2014) show estimates of the near sidelobe 
contribution to the solid angle of between 0.1 and 0.22 %. 

In addition, pre-launch optical calculations show that at a 
few FWHM from the beam centre, the HFI beam will have a 
diffraction pattern from the secondary mirror edges (the pupil 
of the Planck optical system). Figure 14 of Planck Collaboration 
VII (2014) shows that at 353 GHz there are near sidelobes with a 
diffraction signature, which have not been accounted for. These 
features are expected from optical calculations, and they were 
seen in the pre-launch test campaign. A reanalysis of Saturn and 
Jupiter data has been carried out for all frequencies, and the near 
sidelobe wing due to diffraction by the edge of the secondary is 
seen at least at 143 GHz and higher frequencies. Additionally, 
the sidelobe profiles of Planck Collaboration VII (2014) led to 
a removal of the residual time response effects with a weighting 
that also removed optical near sidelobe signal around 20'. 

In order to assess the impact of a more refined accounting 
of the near sidelobes, we have constructed a new hybrid beam 
profile that consists of the nominal beam, azimuthally averaged 
Jupiter observations, and GRASP simulations. The azimuthally 
averaged Jupiter observations extend to 2°, and include all resid- 



Fig. C.2. HFI beam and planet measurement profiles. The grey 
shaded region shows the 1 cr errors in the nominal beam profile 
used to make the effective beam window function. The dashed 
blue and red lines are the best-case and worst-case near side¬ 
lobe estimates from Jupiter, presented in figure 10 of Planck 
Collaboration VII (2014). The solid green curve is the hybrid 
beam profile used to derive correction factors here. The vertical 
dotted lines indicate the radial extent of each portion of the hy¬ 
brid beam: the nominal beam is at small angles; at intermediate 
angles Jupiter data are used; and at large angles the scaled GRASP 
simulation is used. 


ual time response and mirror dimpling effects. The astrophysical 
background is subtracted from each observation, and we fit a 
baseline and a scaling factor to match the Jupiter profile to the 
Mars profile at an overlap angle. The GRASP simulations con¬ 
sist of all eight 100 GHz detectors simulated at five frequencies 
across the band. A 100 GHz full-frequency beam is constructed 
by weighting the simulations with the map-making weights, 
detector bandpass, and planet spectral energy distribution. At 
large angle, the GRASP calculations are dominated by diffrac¬ 
tion around the secondary mirror. The diffraction pattern has a 
similar angular dependence at all HFI frequencies, and its nor¬ 
malization is set by the edge taper of the secondary. We therefore 
scale the 100 GHz simulations to higher frequencies by fitting to 
the Jupiter profile. Figure C.2 shows the radial profiles of the hy¬ 
brid beams, compared with those given in Planck Collaboration 
VII (2014). 

The Saturn-Jupiter combination (Fig. C.2) gives a diffraction 
pattern similar to that predicted and measured preflight (Tauber 
et al. 2010). This increases the estimate of the contribution of 
the near sidelobes to the total solid angle of beam, as given 
in figure 10 of Planck Collaboration VII (2014). However, the 
dominant contribution to the increase in solid angle at 100 and 
143 GHz occurs around 20', and is due to removing the assump¬ 
tion of azimuthal symmetry and applying no time response cor¬ 
rection when computing the radial profile. 

This also affects the shape of the effective beam window 
function, with a decrease similar to that illustrated in the green 
curve in Fig. 19 of Planck Collaboration VII (2014), although 
with a larger amplitude and large error bars (Fig. C.3). The ef¬ 
fect of these corrections for near-sidelobe power missing in the 
beams used in the 2013 results is illustrated in Fig. 12. The error 
in the integral of the beam profile is dominated by the back- 
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Fig. C.3. Effective beam window function corrections based on 
the hybrid beam profile. The green curve shows the correction 
based on the hybrid model presented here. The dashed blue and 
red lines are the best-case and worst-case near sidelobe estimates 
from Jupiter, as shown in Fig. 10 of Planck Collaboration VII 
(2014). The grey shaded region shows the +1 cr errors in the cor¬ 
rection to the solid angle. 


ground removal and inter-calibration of the Mars and Jupiter 
data. 

C.2.2. Far sidelobes 

Estimates from Tauber et al. (2010) put the integrated solid angle 
of the far sidelobes at roughly 0.3 % of the total beam solid an¬ 
gle. While this is small, a significant fraction of this solid angle 
appears near the satellite spin axis. This means that the signal 
seen through this part of the sidelobes is not modulated by the 
satellite rotation, and therefore does not contribute to the calibra¬ 
tion. Using calculations such as those in Sect. C.1.2, we estimate 
the far sidelobe effect on the calibration to be of order one-third 
the total solid angle for the far sidelobes. For HFI, this consti¬ 
tutes about 0.1 % at 100 and 143 GHz. 
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